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intervals of age at exposure and time after exposure were used in these
analyses. Because the numbers of cases of cancer were relatively small, and
because the risk of colon cancer may be related to ankylosing spondylitis
itself, analyses were restricted to leukemia and, as a group, all other cancers
except colon cancer.

Statistical Methods

The program AMFIT, described in Annex 4C, was used to fit various
exposure-time-response models to these data sets. This program fits a
general form of "Poisson regression" model, in which the observed number
of events in each cell of the cross-tabulation is treated as a Poisson variate
with parameters given by the predicted number of events under the model,
the product of the person-years in that cell times the fitted rate. The specific
models used can be formally expressed as follows. Let 70 denote the age-
specific background risk of death due to a specific cancer for an individual
at a given age. This background risk will also depend upon the individual's
sex and birth cohort (that is year of birth). For a given radiation dose
equivalent d in sievert (Sv) we write the individual's age-specific cancer risk
j(d) as

7(d) = 7o[l + /(<W)].                           (4-1)

Let/(rf) represent a function of the dose d which in the committee's models
is always a linear or linear-quadratic function, i.e., f(d) = aid or f(d) =
a^d + <*3<i2. In general, the excess risk function, g(/3) will depend upon a
number of parameters, for example, sex, attained age, age-at-exposure, and
time-since-exposure. One can also write the age-specific risk as an additive
risk model

7(rf) = 7o + /(<W)-                              (4-2)

These models give similar results (see Annex 4D) as expected since the
function g(/3) is allowed to depend on age, time, etc. This would not be
the case if g(/?) were restricted to having a constant value other than for
sex and age at exposure.

The models were fitted using maximum likelihood, i.e., the values of
the unknown parameters which maximize the probability of the observed
number of cases (the "likelihood function") are taken as the best estimates,
and, where applicable, confidence limits and significance tests are derived
from standard large-sample statistical theory.

It was expected that the form of the background term might vary
considerably between populations at risk and is not of particular interest in
terms of radiation risk. The committee chose not to model it, but rather